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ABSTRACT 

"Extreme Scattering Events" (ESEs) are attributed to radio-wave refraction by a 
cloud of free-electrons crossing the line-of-sight. We present a new model in which these 
electrons form the photo- ionized 'skin' of an underlying cool, self-gravitating cloud in the 
Galactic halo. In this way we avoid the severe over-pressure problem which afflicts other 
models. The UV flux in the Galactic halo naturally generates electron densities of the 
right order. We demonstrate, for the first time, a good reproduction of the prototypical 
ESE in the quasar 0954+658. The neutral clouds are a few AU in radius and have masses 
< 1CT 3 M . The observed rate of ESEs implies that a large fraction of the mass of the 
Galaxy is in this form. 

Subject headings: scattering — ISM: clouds — Galaxy: halo 
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1. Introduction 

ESEs were discovered ten years ago during radio 
flux monitoring of a sample of compact radio quasars 
(Fiedler ct al 1987). The ESE phenomenon consists of 
dramatic flux changes occuring over several weeks to 
months. It is broadly agreed that ESEs are not intrin- 
sic variations but, rather, apparent flux changes which 
are caused by refracting elements, a few AU in ra- 
dius, crossing the line-of-sight. Both random (Fiedler 
et al 1987) and deterministic (Romani, Blandford & 
Cordes 1987) lens structures have been proposed; in 
all cases the refraction is attributed to free electrons. 
Two further points of consensus are that the blobs 
of free electrons must be Galactic, and that they 
represent a distinct component of the ISM (Narayan 
1988; see also Rickett 1990). This same component 
is thought to be responsible for episodes of multiple 
imaging of radio pulsars (Cordes & Wolszczan 1986; 
Rickett 1990) — a phenomenon which manifests itself 
as periodic fringes in the dynamic spectra. 

Despite the attention which the ESE phenomenon 
has attracted (Romani, Blandford & Cordes 1987; Ro- 
mani 1988; Clegg, Chcrnoff & Cordes 1988; Clegg, 
Fey & Lazio 1998), the current state of understand- 
ing is unsatisfactory for two reasons. First the im- 
plied pressure of the electron cloud - having an in- 
ferred density n e ~ 10 3 cm~ 3 and a temperature 
T ~ 10 4 K - exceeds that of the diffuse Interstellar 
Medium (ISM) by a factor of 10 3 . Such a cloud should 
dissipate on a time-scale comparable to the sound- 
crossing time, of order a year. This difficulty can be 
lessened by appealing to cloud geometries which are 
elongated along the line-of-sight, and by placing the 
clouds in regions of high pressure within the ISM — 
specifically, old supernova shocks. In support of the 
latter idea it has been argued (Romani 1988) that the 
distribution on the sky of recorded ESEs reflects that 
of the strongest regions of non-thermal radio emission. 
However, with the small number of ESEs recorded to 
date, and the large area covered by the radio loops, 
this evidence is at best suggestive. The second main 
difficulty with existing models is that none has been 
able to reproduce the dual-frequency light-curve of 
the quasar 0954+658 — the source for which the best 
ESE data exist. 

In this paper we present the essential elements of a 
new model for ESEs, wherein the free electrons com- 
prise a 'skin' of photo-ionized material around a cool, 
self-gravitating cloud. The pressure of the neutral 



material is balanced by the self-gravity of the cloud, 
while the ionized gas flows continuously from the sur- 
face; complete evaporation occurs over a time-scale 
of order the Hubble time. We show in §2 that this 
model reproduces the prototypical ESE light-curves. 
The event rate for ESEs then leads to the robust con- 
clusion (§3) that these clouds contribute substantially 
to the Galactic dark matter. In §4 we note some key 
observational tests. 

2. Refraction in a photo-evaporated wind 

The estimated cloud dimension and distance indi- 
cate that geometric optics is a good approximation. 
To proceed we need to know the electron density dis- 
tribution, from which we can calculate the bending 
angle of rays, and hence the imaging properties of the 
cloud. 

The photo-ionized material at the surface of the 
cloud is not bound, but flows away as a feeble wind. 
Dyson (1968) constructed an analytic solution for 
the density profile of a spherically symmetric photo- 
evaporated wind from a neutral cloud; his solution is 
plotted in figure 1. The electron density at the surface 
of the cloud is given (in cm -3 ) by 4.87 x 10 6 y/Joo/R, 
and in the local Galactic halo, at the solar circle, the 
ionising photon field is Joo — 10^ cm 2 s * sr * (Dove 
and Shull 1994). For a cloud radius fi~3x 10 13 cm 
(Fiedler et al 1987) we obtain n e (R) ~ 10 3 cm" 3 . 
This model thus exhibits immediate strength in that 
the electron density required to explain the ESEs is a 
prediction of our model, whereas it is an assumption 
of existing models. 

Dyson's solution cannot be used directly for our 
calculations, however, as it is discontinuous at the 
ionisation front, whereas we require the first and sec- 
ond derivatives of n e to be well behaved. Physically 
the thickness of the ionisation front is just the ab- 
sorption length of Lyman-limit radiation in the neu- 
tral medium. In turn this is fixed by the tempera- 
ture of the neutral atmosphere, via the requirement 
for pressure balance across the transition region. An- 
other difficulty with Dyson's solution is that n e is 
determined implicitly, which is computationally in- 
convenient. We therefore employ the approximation 
n e (x) oc F(x) x~ 2 (l + 41og e re) -1 / 2 , where x is the ra- 
dius expressed in units of the cloud radius, and F is 
the Fermi function, taken to have a width of fa = 0.01 
(corresponding to a temperature of 40 K in the neu- 
tral atmosphere). This distribution is plotted in figure 
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1. 

Making use of this electron density distribution, 
together with the phase-screen approximation, it is 
straightforward to compute light-curves for a neutral 
cloud traversing the line-of-sight to a distant radio 
source. Figure 2 shows one such computation, appro- 
priate to a cloud at a distance of 2 kpc and an impact 
parameter of R/2, along with the data for 0954+658. 
To simulate the non-zero size of the real source, the 
theoretical curves are smoothed by convolution with 
Gaussian functions. The adopted source model cor- 
responds to 50% of the source flux being contained 
within a compact (lensed) component, with this com- 
ponent having a brightness temperature of 8 x 10 11 K. 

The main qualitative features of these light curves 
are accounted for as follows. The phase velocity of 
the wave is increased by the presence of free elec- 
trons, so the cloud acts as a diverging lens. At low 
frequencies the lens is powerful enough that almost 
all rays are refracted out of the line-of-sight, and only 
a small flux is measured when the lens is aligned with 
the source; this behavior is generic to all blobs of free 
electrons regardless of the details of their density dis- 
tribution. Consequently this regime of a very strong 
lens is not particularly helpful in distinguishing our 
model from other possible electron density distribu- 
tions. At higher frequencies, however, the refractive 
index of the lens is much smaller, and typically rays 
are no-longer refracted through sufficiently large an- 
gles that caustics form. Exceptions occur for rays 
that pass near the edge of the cloud: here the pho- 
toionized skin creates a large phase curvature in the 
wavefront and can introduce caustics even at high fre- 
quencies. These caustics are evident as sharp peaks 
in the model light-curve, and we note that similar 
features appear in the data; in our model they are 
intimately associated with the presence of a peak in 
the electron column density at the limb of the cloud. 
The good qualitative agreement between our model 
and the data suggests that surface photo-evaporation 

- which, of course, implies underlying cool material 

- is an essential feature of real ESE clouds. By con- 
trast, the Gaussian electron density profile originally 
proposed (Romani, Blandford and Cordes 1987) for 
ESEs cannot, even qualitatively, match the dual fre- 
quency light curve of 0954+658. 



3. Implications for dark matter 

Going beyond the interpretation of individual events, 
the principal implication of the new model is that 
there is much more mass present in the ESE clouds 
than had been previously thought (cf. Fiedler et al 
1987; but see also Pfenniger & Combes 1994). We 
now derive a lower limit on the total contribution of 
the ESE clouds to the mass of the Galaxy. 

The sky-covering factor, /, of the clouds is esti- 
mated from the flux monitoring data (Fiedler et al 
1994) to be / — 5 x 10~ 3 . We can immediately re- 
late / to the total surface density, S, in clouds at the 
solar circle: £ ~ 2 (sin |b|)/ M/irR 2 ~ / M/irR 2 , for 
cloud mass M and radius R (ESEs do not preferen- 
tially occur at low Galactic latitude, b, so we have 
set (sin |b|) ~ 0.5). These quantities arc in turn re- 
lated by the requirement of hydrostatic equilibrium 
within each cloud - fcT ~ GMm p / R at temperature 
T - leading to £ ~ fkT/TrGm p R. The cloud radius 
can be inferred from the event duration, in combina- 
tion with an assumed transverse speed; a limit fol- 
lows from setting the transverse speed equal to the 
escape speed for the Galaxy (500 kins" 1 ), giving 
R < 3 x 10 14 cm. We expect that the cloud tem- 
peratures are at least as large as that of the cosmic 
microwave background (i.e. T > 3 K), so we deduce 

S > 10 2 Mq pc- 2 . (1) 

This lower limit is already larger than the total sur- 
face density of observable matter, and consistent with 
the value (210 M Q pc~ 2 ) necessary to explain the 
Galactic rotation curve (Binney and Tremaine 1987). 
Taking the dynamically determined surface density 
as an upper limit on the contribution from ESE 
clouds, we estimate the individual cloud masses to be 
< 10~ 3 M Q : highly sub-stellar. Hence the ESE data 
lead us to the conclusion that low-mass gas clouds 
make up a good fraction of the mass of the Galaxy. 

We are currently reassessing the sky covering frac- 
tion, / (manuscript in preparation) , based on our the- 
oretical model and the data of Fiedler et al. (1994). 
It is possible that / is as small as a few times 10~ 4 , in 
which case the lower limit on £ would be an order of 
magnitude smaller. Even at this level, however, the 
cool, low-mass cloud population implied by the model 
remains a substantial and poorly understood compo- 
nent of the Galaxy. Moreover, because our estimate is 
a lower limit on the surface density, a smaller covering 
fraction still admits the possibility that such clouds 
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do indeed dominate the dynamics of the Galaxy. 

4. Discussion 

The result we have presented might appear sur- 
prising; can dark matter really be in this form and 
yet have escaped detection? It has been noted pre- 
viously (e.g. Pfenniger, Combes & Martinet 1994, 
Gerhard & Silk 1996) that cool, dense gas clouds are 
not easy to detect directly and may indeed be a vi- 
able dark matter candidate. We reconsider these is- 
sues elsewhere (manuscript in preparation), but some 
brief comments are appropriate here. 

First, by virtue of their high density (~ 10 12 cm~ 3 ) 
and low temperature, three body reactions are suffi- 
cient to ensure that the clouds are molecular. The 
constituent material must have a very low dust/gas 
ratio in comparison with the diffuse ISM, otherwise 
the population would already have been revealed 
by optical extinction events of extragalactic stars 
(notably in the data generated by the microlensing 
searches: Paczyhski 1996). A low dust content might 
reflect either a very low metal abundance or, possi- 
bly (B. Paczyhski, personal communication), settling 
of the refractory elements into a small, rocky core. 
Thirdly, unless the clouds are dynamically unimpor- 
tant (i.e. S <C 10 2 M Q pc~ 2 ), the measured 7-ray 
background at high Galactic latitudes (Kniffen et al. 
1996) constrains the vertical distribution of the clouds 
to be much more extensive than the cosmic ray disk, 
implying that the clouds are a halo population. (In 
turn this means that they are pre-Galactic objects.) 
The small covering fraction, /, of the clouds means 
that each can make 1// orbits of the Galaxy before 
colliding with another cloud, implying a lifetime of 
order the Hubble time. The large column density 
of the individual clouds (of order 10 25 cm -2 ) means 
that interactions with the diffuse gas in the Galac- 
tic plane involve a negligible deceleration. The life- 
time against UV photo-evaporation is also of order 
the Hubble time. 

While the cool-cloud model we advance for ESEs 
is compatible with a variety of existing data, some 
specific observations can be made which would pro- 
vide strong tests of the picture. We offer three such 
tests; the first of these is most intimately connected 
with ESEs, and is in some sense the easiest experi- 
ment to perform. Although the clouds are predomi- 
nantly molecular, a crude assessment of the relevant 
chemistry suggests that they should contain sufficient 



atomic hydrogen to make them opaque in the 21 cm 
line. Thus, while an ESE is in progress there should 
also be a narrow absorption line present, in addition 
to the normal absorption spectrum of the source. (In 
our model, a similar connection should also hold be- 
tween 21 cm absorption and multiple imaging of pul- 
sars.) Small-scale (tens of AU) structure has already 
been reported in 21 cm absorption lines (Davis, Dia- 
mond & Goss 1996, Frail et al 1994), but the possi- 
ble connection with ESEs has not yet been explored. 
The line velocities should reflect the kinematics of the 
cloud population, which are expected to be quite dif- 
ferent to that of the diffuse HI in the Galactic disk. 
We note that, because of the multiple imaging which 
takes place, the absorption lines may appear satu- 
rated but with a non-zero flux minimum. 

Secondly, because the inferred column density of 
neutral material is so large it is expected that the 
clouds are opaque to Rayleigh scattering (c.f. Combes 
& Pfenniger 1997) between A ~ 2000 A and the Ly- 
man/ Werner bands of H 2 ; and opaque to Thomson 
scattering at X-ray energies. Thus one in 1// ex- 
tragalactic UV/X-ray sources (of sub-milliarcsecond 
size) should appear extinguished at any one time, 
with events lasting a month or so. 

Finally we note that the recombination lines (Ha 
etc) from the cloud surface may be detectable if there 
is a sufficiently strong local source of ionising radia- 
tion (Gerhard & Silk 1996). That is, clouds which are 
within a parsec of a hot star will be luminous sources 
of Ha and may thus be rendered visible. It might 
also be possible to detect, from their Ha emission, 
the clouds which are closest to the Sun. Many clouds 
are expected within a parsec, and these may appear 
as compact Ha nebulae of large radial velocity and 
very high proper-motion. 

5. Conclusions 

We have presented a new model for the cause of 
Extreme Scattering Events, where radio-wave refrac- 
tion by ionized material magnifies ( "lenses" ) a distant 
radio source. In our model the ionized material is 
generated by photo-evaporation of an underlying neu- 
tral, self-gravitating cloud, so there is no difficulty in 
understanding the requisite high electron pressures; 
indeed the ambient UV field in the Galactic halo nat- 
urally generates the necessary electron density. More- 
over we find it straightforward to reproduce the dual- 
frequency light-curve of the prototypical ESE in the 
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source 0954+658. These facts argue strongly for the 
soundness of the model. 

It follows immediately that a large fraction of the 
mass of the Galaxy is in the form of cool gas clouds of 
up to Jovian mass, and only a few AU in radius. We 
are unable to falsify this picture with existing data, 
but with new observations some straightforward tests 
can be made. 

Particular thanks go to Ron Ekers for many helpful 
discussions; we are also grateful to Simon Johnston, 
Barbel Koribalski, Don Melrose, Bohdan Paczynski, 
Gordon Robertson and Lister Staveley-Smith for their 
input. 
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Fig. 1. — Electron density as a function of radius 
for a photo-evaporated wind: Dyson's (1968) solution 
(dashed line); and our approximation (full line). 



This 2-column preprint was prepared with the AAS IATj^X 
macros v4.0. 
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Fig. 2. — Observations (top panel, adapted from 
Clegg, Fey and Lazio [1998]; original data in Fiedler 
et al. [1987]) and theory (lower panel) for the ESE 
seen in the quasar 0954+658. The upper curves corre- 
spond to a radio frequency of 8.1 GHz (to which 1 Jy 
has been added), and the lower curves to 2.7 GHz. 
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